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To expand the framework available for interpreting experiments on disordered strongly correlated 
systems, and in particular to explore further the strong-coupling zero-bias anomaly found in the 
Anderson-Hubbard model, we ask how this anomaly responds to the addition of nonlocal electron- 
electron interactions. We use exact diagonalization to calculate the single-particle density of states 
of the extended Anderson-Hubbard model. We find that for weak nonlocal interactions the form of 
the zero-bias anomaly is qualitatively unchanged. The energy scale of the anomaly continues to be 
set by an effective hopping amplitude renormalized by the nonlocal interaction. At larger values of 
the nonlocal interaction strength, however, hopping ceases to be a relevant energy scale and higher 
energy features associated with charge correlations dominate the density of states. 
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I. INTRODUCTION 



As the simplest model which incorporates strong cor- 
relations, the Hubbard model is widely used as a start- 
ing point for exploring the diverse behaviors of transi- 
tion metal oxides. Understanding the effect of disor- 
der in these systems is of interest given the importance 
of chemical doping in tuning their properties. In fact, 
nanoscale electronic disorder has been observed and ex- 
plicitly correlated with the locations of dopant atoms in 
some materialsii Disorder can be incorporated into the 
Hubbard model, and the most widely-studied version is 
the Anderson-Hubbard model (AHM), which is defined 
below. The combination of local interactions and disor- 
der in the Anderson-Hubbard model has been shown to 
result in a suppression of the density of states (DOS) at 
the Fermi level, a feature known as a zero-bias anomaly 
(ZBA). The ZBA in the AHM has a number of unique 
features: The DOS suppression occurs in the absence of 
nonlocal interactions, in contrast with other, well-known 
examples including both the Efros-Shklovskii Coulomb 
gap^ and the Altshuler-Aharonov anomalyiS.. Moreover, 
the energy scale in the AHM ZBA is proportional to the 
hopping amplitude t and independent of both the local 
interaction U and the disorder strength A within a broad 
range of phase space. 

The strong-coupling ZBA appears to arise from kinetic 
energy savings, rather than a reduction of the Coulomb 
energyj^ In the absence of hopping, electronic states are 
confined to individual atomic sites. With nonzero t, these 
states can extend over nearby sites, lowering the elec- 
tronic kinetic energy. In the large U limit, however, this 
spreading of electronic wave functions is strongly inhib- 
ited by the local Coulomb interaction. Nonetheless, if 
the disorder is sufficiently strong, there will be atomic 
sites for which the energy of double occupancy is nearly 
degenerate with the energy to form a singlet with one of 



its neighbors. For these special configurations, the local 
Coulomb interaction does not inhibit electronic motion 
between the two atoms, and the energy of the system is 
reduced by an amount of order t relative to the atomic 
case.i^ This leads directly to the suppression of spectral 
weight at the Fermi energy over an energy scale t. This 
mechanism is unique to strongly-correlated systems. 

Given the presence of nonlocal interactions in real ma- 
terials, as well as their importance in theories of ZBAs 
in other models, it is natural to ask how nonlocal in- 
teractions influence the kinetic-energy driven ZBA in 
the AHM. This formal question bears on a number of 
classes of materials. (1) Most directly this work relates 
to doped transition metal oxides. DOS measurements^ 
in SrRui-xTia^Os and LaNii-xMn^Oa show deviations 
from the standard pictures of Efros-Shklovskii^ and 
Altshuler-Aharono\«^, and it is an open question whether 
this is because of strong correlation physics. In addi- 
tion, early work on the ZBA in the AHM suggested that 
disorder may contribute to the stability of the pseudo- 
gap in high temperature superconductors!^ (2) A second 
class of materials to which this work may be relevant are 
dilute doped semiconductors and granular metals. Gen- 
erally described using atomic-limit models, these systems 
display Coulomb gap behavior. Open questions include 
whether there is an association between the Coulomb gap 
and glassy behavior;^"— and how electron mobility (and 
hence screening) influence the observed behaviors »iiii^ 
(3) A third class of materials are two-dimensional elec- 
tron gas systems, such as thin metal films and MOS- 
FETs. Whereas (noninteracting) localization theory con- 
cludes that there should be no metal-insulator transition 
in two dimensions, experiments on these systems sug- 
gest otherwiseJ^ Recently it has been proposed that the 
insulating behavior in these films is in fact due not to 
disorder but to interactions.— To support this proposal, 
the extended Hubbard model was explored, but only in 
the clean limit. Even if interactions drive the insulating 



behavior, the systems remain disordered. (4) Finally, or- 
ganic conductors are another class of materials in which 
strong correlations can be important. Very recently it 
has been shown that disorder can be introduced into or- 
ganic conductors by x-ray irradiation, resulting in novel 
behaviors and expanded opportunities for exploring in- 
teractions and disorder ji^ii^ 

To expand the theoretical framework available for in- 
terpreting these diverse materials, we explore the DOS 
of the extended Anderson- Hubbard model (EAHM) on a 
number of trajectories in the available phase space, all 
in the strong disorder limit. Figure [T] summarizes our 
main result, namely the presence of a crossover in the 
form of the DOS as a function of interaction strength. 
When the nonlocal interactions are weak, electron mobil- 
ity plays a key role, generating in particular the narrow 
kinetic-energy-driven ZBA seen earlier in the AHM.— i^ 
In the half-filled case shown in Fig. [I] when nonlocal 
interactions are strong, atomic-limit physics dominates 
the DOS: A broad suppression around the Fermi level is 
associated with charge correlations driven by the nearest- 
neighbor repulsion. Doping away from half filling reduces 
the impact of interactions and hence attenuates this ef- 
fect. 

Section In] describes our approach while our results are 
presented and discussed in Section IlIII 
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FIG. 1: Density of states versus frequency at half filling {n = 
1) for A = 12i, U — 8t and V values as indicated. Throughout 
this work, results are for 12-site lattices averaged over 1000 
disorder configurations. Inset indicates the location of each 
parameter set on the phase diagram of the clean extended 
Anderson-Hubbard model. 



II. METHOD 

The extended Anderson-Hubbard model includes the 
hopping t and onsite Coulomb repulsion U of the tradi- 
tional Hubbard model as well as two additional terms: 
a nearest-neighbor Coulomb repulsion V and disordered 



site potentials e^. 

y- = -t Yl Ccj<T + Yl f^^'t"4 

(ij>-cr i 



(1) 



We consider a two-dimensional square lattice, c]^ is the 
creation operator for lattice site i and spin a. fiia = 
cj^Cia- {i,j) refers to nearest neighbor pairs. The site 
potentials e^ are chosen from a flat distribution of width 
A: P(e,) = e(A/2 - |ej|)/A where 9 is the Heaviside 
function. We focus on the limit of strong disorder: All 
results shown are for A = 12 in units of the hopping t. 

The Lanczos method is used to calculate the DOS of 
12-site clusters. The Lanczos methodic denotes a collec- 
tion of iterative procedures all founded on the idea that 
a matrix Q can be found such that Q iL Q — 21 where 
T is a tridiagonal matrix. Computational savings come 
from the fact that the number of columns n in Q may 

be less than the number of rows such that T is smaller 
than _ff . The extremal eigenvalues of T converge quickly 
to those of _ff as a function of increasing n. Finding the 
DOS using the Lanczos method proceeds in two steps. 
First, restarted Lanczos is used to find the ground state 
I'I'o) and energy Eq- Second block- recursior>i^ is used to 
calculate the Green's function. 

The LDOS at site i of a particular disorder configura- 
tion c is given by 
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is the r diagonal element of the real-space Green's func- 
tion. Here |V'o) and Eq are the ground state wave func- 
tion and the ground state energy of disorder configuration 
c. The DOS of a single disorder configuration is 
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where A^^ = 12 is the number of sites in the lattice. We 
present DOS results averaged over many disorder config- 
urations 
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The number of disorder configurations Nconfig ~ 1000 
for all results presented here. Because we study systems 
with very strong disorder, the mean free path is of order 
the lattice spacing, and hence the disorder averaged DOS 
can be expected to be representative of the DOS in the 
thermodynamic limit. Comparison of Lanczos results on 
a 10-site cluster with determinant quantum Monte Carlo 
results on a 64-site system by Chiesa, et al^ support this. 



III. RESULTS AND DISCUSSION 

To explore how nonlocal interactions affect the ZBA 
in the Anderson-Hubbard model, we present here DOS 
results first for the case of half filling, followed by lower 
dopings. 



A. Half filling 
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FIG. 2: DOS versus frequency at half filling (n = 1) for 
A = 12 and [7 = 8. (a) Analytic result for 1/ = and 
4 = with the occupancies of contributing sites indicated, 
(b) Diagram showing an arrangement of atomic orbitals for 
which the introduction of hopping moves DOS weight away 
from the Fermi level, (c) Numerical results for V = and f 
values as indicated, (d) Numerical results for Y = 1.6 and f 
values as indicated. 



1. Review of the V = case 

Before proceeding to the new results, let us review 
what is known about the V = Q case. The highest curve 
in Fig. [1] shows the ZBA found in the AHM, which was 
explored by Chiesa, et al.— For sufficiently strong disor- 
der, the width of this ZBA is independent of disorder 
strength, interaction strength, and doping. The one pa- 
rameter which controls the width is the hopping i, on 
which the width depends linearly. The existence of a 
ZBA in the absence of nonlocal interactions is novel. The 
two standard frameworks for discussing ZBAs were de- 
veloped by Efros and Shklovskii^ and by Altshuler and 
Aronov^ Efros and Shklovskii^ addressed a system with 
l/r Coulomb interactions and disordered site potentials 
in the atomic limit. They argued that the DOS at an 
energy e is proportional to |e — efI''""^ where ep is the 
Fermi level and d is the dimension. The zero DOS at 
the Fermi level is known as a soft gap or more specifi- 
cally as the Coulomb gap, and its existence depends on 



the infinite range of the Coulomb interaction. Altshuler 
and Arono\*^ considered the limit of weak disorder and 
weak interactions, and used diagrammatic perturbation 
theory to show that in this limit a cusp appears in the 
DOS near the Fermi level. The result is not especially 
sensitive to the form of the interaction. However, for 
purely local interactions the correction to the DOS is 
positive; DOS suppression only occurs when nonlocal in- 
teractions are present. In both the Altshuler-Aronov and 
Efros-Shklovskii pictures nonlocal interactions are key to 
the suppression of the DOS at the Fermi level, yet in the 
AHM there is a ZBA. Moreover, although a dependence 
on the effective nonlocal interaction J oc t^ /U might be 
expected, here the dominant energy scale is t alone. 

To understand this linear dependence on f, it is useful 
to start from the atomic limit and then consider what 
happens as hopping is turned oni^i^ '^^i^" In the atomic 
limit, each site contributes to the DOS at, at most, two 
energies: the site potential e^, and et + U. It is convenient 
to refer to these as the lower Hubbard orbital and the up- 
per Hubbard orbital. When the site is singly occupied, it 
contributes to the DOS at both energies because a par- 
ticle may be either added or removed. When the site is 
empty (doubly occupied) only the addition (removal) of 
a particle is possible and hence the DOS contribution is 
only at Ci (e^ + U). In the ground state, sites with poten- 
tials ei > ^ are empty, those for which fi > ei > ^ — U are 
singly occupied, and those for which ^— U > ei are dou- 
bly occupied. Putting these together, the atomic-limit 
DOS corresponding to any combination of A, U and fi 
may be constructed. The case of A = 12, t/ = 8 and 
/i = 4 (half filling) is shown in Fig. [2Ja). The numbers 
in each block indicate the ground state occupancy of the 
sites which contributed. An important point for our pur- 
poses is that there is no ZBA at zero temperature in the 
atomic limit. Nonzero temperature docs suppress the 
DOS even in the atomic limitr^ but we restrict ourselves 
here to zero temperature. 

We now ask how hopping affects the DOS. To address 
this it is convenient to consider the simple case of a 
two-site system. Consider in particular the configura- 
tion shown in Fig. ^h): ei + U is just below the chemical 
potential, and £2 is just above. In the atomic limit, the 
first site will be doubly occupied and the second empty, 
corresponding to the Fock state |20) and to a grand po- 
tential E — jiN near zero. However, when hopping is 
allowed, the new ground state will be a linear combina- 
tion of the singlet states (| tl) - I lt))/v^ and |20). The 
grand potential of this ground state is lower than the 
atomic one by t, because the probability amplitude for 
finding an electron is now spread over both sites lowering 
the kinetic energy of the many-body state. To linear or- 
der in f, there is no corresponding shift in the 1-particle 
and 3-particle excited states. The energy of transitions is 
therefore increased, corresponding to a shift of the poles 
in the Green's function away from the Fermi level. If we 
consider an ensemble of such two-site systems, the result 
is the opening of a ZBA in the DOS of width t^^ 



While larger lattices present many additional 
complications;^ this simple two-site picture demon- 
strates how kinetic-energy savings can lead to a ZBA of 
width t. It is worth emphasizing that this behavior is 
unique to strongly correlated systems because it relies 
on there being a large difference in energy between the 
atomic lower and upper Hubbard orbitals at each site. 
For this reason, the effect is not captured by mean-field 
treatments. Finally, we note that the effect requires 
double occupancy on some nonzero fraction of sites. Fig. 
[5Ja) shows that, in the atomic limit when both empty 
and doubly occupied sites are present, the Fermi level 
falls somewhere within the high central-plateau. When 
the chemical potential is lowered such that no sites are 
doubly occupied, the Fermi level instead sits at the left 
edge of this central plateau. When hopping is nonzero 
but still much less than U and A. this feature of the cen- 
tral plateau in the DOS persists, and the position of the 
Fermi level relative to this central plateau continues to 
be indicative of the level of double occupancy. Because 
of the importance of double occupancy in the formation 
of the ZBA, the linear t dependence of the ZBA is only 
expected when the Fermi level falls within this central 
plateaui^ 



2. Evolution of the DOS with V 



at t = is a manifestation of the same atomic limit 
physics found in the Efros-Shklovskii Coulomb gap, but 
for a short-range interaction. The DOS is not suppressed 
to zero here because the interaction range is finite. 

When hopping is turned on, the shape of the ZBA 
changes abruptly. Note that the t ^ 0.5 curve coincides 
with the t = 1 curve at low energies and with the t ~ 
curve at higher energies. This suggests that the reshap- 
ing of the ZBA by hopping begins at the Fermi level and 
spreads outward in energy as t is increased. The atomic 
limit (t = 0) is classical in the sense that only integer 
occupancy is allowed. It appears that the quantum ef- 
fects introduced by hopping have their first effect on the 
DOS at the Fermi level, while classical behavior persists 
at higher energies. 

In summary. Fig. [5] emphasizes the sharp distinction 
in the atomic limit between having nonlocal interaction 
and not. The addition of hopping generates a new ZBA 
which is qualitatively similar with and without nonlocal 
interactions. 



4. Small V 



Having reviewed the V = case, we return to Fig. [T] 
The remaining curves demonstrate the evolution of the 
DOS as a function of the nearest-neighbor interaction 
strength V. Qualitatively, the lower-]/ curves resemble 
the V ~ curve and arc distinct from those at higher 

V values. For V = 0.8 and V ~ 1.6 the most obvious 
changes are in the height and width of the central peak, 
while the form of the ZBA is relatively consistent. For 

V = 2A and V = 3.2, however, there is a strong shift 
of spectral weight away from the Fermi level, such that 
the DOS is largest near the band edges. Moreover, the 
ZBA loses its sharp form and opens into a hard gap. The 
sections below provide a more detailed discussion first of 
the distinct physics present when y 7^ 0, then of the 
small V behavior and finally of the large V behavior. 



3. Comparing V = with V ^ 

A striking feature of Fig. [1] as noted above, is the 
similarity in the shape of the ZBAs seen at T^ = 0.8 and 
1.6 with that at ]/ = 0. Fig. [2] emphasizes an important 
distinction between the case of V^ = and that oiVy^O. 

For the case V = 0, Fig. ^c) shows the DOS with and 
without hopping. The t = curve shows the wedding- 
cake structure predicted in Fig.[5Ja), with no ZBA. When 
t is turned on, a ZBA emerges with a width linear in t as 
described above. 

In contrast, when V ^ 0, Fig. ^d) shows that there is 
a ZBA even in the atomic limit. This rounded anomaly 




FIG. 3: Symmetrized DOS versus frequency in units of t at 
half filling (a) and at quarter filling (c) shifted to coincide 
at the Fermi level. A — 12t, U = 8i, and V values as indi- 
cated. Panels (b) and (d) show the same data plotted versus 
frequency in units of the mean-field rescaled hopping i. 



We now turn to a more quantitative comparison of the 
three qualitatively similar curves: ^ = 0, 0.8 and 1.6. 
Fig. [3] (a) shows the three curves vertically shifted to co- 
incide at the Fermi level. Here we see that the ZBA is 
becoming broader as V is increased. For these small V 
values, a mean- field treatment of the nonlocal interac- 
tion provides some insight. A mean-field treatment of 
the nearest-neighbor interaction results in the following 



expression: 

V 
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where Ui = J^ai^a^ia) and fjia = (Cj^Ci^) with (...) de- 
noting the expectation value with respect to the ground 
state. The last two terms are constants which simply 
shift the zero of energy. The first term in this expression 
renormalizes the site energies e^ consistent with the idea 
that interactions screen the disorder potential. However, 
when V is very small relative to the disorder potential 
and U is large such that the charge density is very uni- 
form, this term will have very little effect. Our focus 
here is on the second term which results in a renormal- 
ized hopping integral: 

where i and j are nearest neighbor sites. This is con- 
sistent with the gradual increase in the width of the 
anomaly as V is increased. 

More precisely, Fig. [3l^b) shows the V = 0, 0.8 and 1.6 
curves with the frequency axis in units of t. Under this 
rescaling, the V^ = 0.8 curve coincides very closely with 
the V^ = curve in the frequency range of the ZBA. The 
V^ = 1.6 curve also matches but only at very low energies. 
That this mean- field approach works less well at V = 1.6 
than at V = 0.8 is due in part to the other changes V 
causes in the system as discussed further in the next sec- 
tion. However, another issue here is that in calculating 
t for this figure we used /,,,■ averaged over all bonds in 
the lattice. With more computational effort, a value for 
t more specifically associated with the sites which con- 
tribute to the DOS near w = could be constructed. We 
expect fij to be larger on bonds between sites with or- 
bital configurations as in Fig. [Ub) , and hence we expect 
that this improved mean-field treatment would create a 
stronger rescaling resulting in a better match at V ~ 1.6. 



5. Large V 

Interaction strength generally falls off with distance, 
making V values near U unphysical. However, well before 
this cutoff, qualitatively distinct behavior arises: The 
V — 2A and 3.2 curves in Fig. [I] are very different from 
those at lower V values. There is a strong shift of spec- 
tral weight well away from the Fermi level such that the 
DOS is largest near the band edges. 

A useful point of reference is the clean extended Hub- 
bard model. In two dimensions, this model has a 
first-order phase transition between charge density wave 
(CDW) order for U < W and spin density wave (SDW) 
order for U > AV^ Fig. [T] inset shows this clean phase 



diagram with the locations corresponding to the DOS 
curves marked. 

How does disorder affect this phase diagram? 
Quenched disorder influences first-order phase transitions 
in a wide variety of ways depending on the details of the 
model.— One possibility is that the phase boundary may 
be moved. However, in our case the boundary is between 
two ordered phases, neither of which is enhanced by the 
disorder. Indeed, the fact that the abrupt change in the 
shape of our DOS occurs between T^ = 1.6 and V = 2.4 
suggests that U = AV remains significant. The bigger 
issue is whether regions of order remain on either side of 
this line. 

On the [/ = axis, a connection may be made 
with the well studied random-field Ising model (RFIM). 
The atomic limit of the EAHM at f/ = is "H = 
J2i''t'^t + ^J2{i.j)'>T-t'^j- Using Hi = 5.1 + 1, this be- 
comes v. = J2i ^iSi + VJ2{i j) SiSj+ constants (for fixed 
particle number). When U — 0^ site occupancies rii are 
restricted to and 2 at zero temperature. Therefore, 
this is precisely the RFIM with ci playing the role of the 
local field and V the spin interaction favoring antiferro- 
magnetic (AFM) order. In two dimensions, the RFIM 
is always disordered!^ Although a more rigorous proof 
has been developed^!, this is most easily illustrated by 
the following surface-to- volume argument i^ Imagine the 
system begins with perfect AFM order, and then con- 
sider flipping all the spins in a domain of size L'^ where 
d is the dimension. Such a flip will raise the energy as- 
sociated with the interaction term by an amount propor- 
tional to the length of the boundary: L'^~^. This flip will 
also change the energy associated with the random field. 
By the central limit theorem, this change in field energy 
has an average value of zero and root mean square value 
proportional to L'^^^. When the field energy savings are 
less than the interaction energy cost (L''/^ < L''-^^) the 
system remains ordered. However, for d < 2, domain 
formation is favored, and the system is disordered for 
any nonzero disorder strength. Smaller ratios of disor- 
der strength to interaction strength correspond to larger 
characteristic domain sizes. 

Returning to the EAHM, the correspondence with the 
RFIM shows that in the atomic limit and with f7 = the 
EAHM will have no CDW order for any nonzero disor- 
der strength. Moreover, hopping and onsite interactions 
both lower the energy cost of the boundaries. Hopping 
makes the site occupancies continuous variables which 
can vary smoothly across domain walls. Onsite inter- 
action promotes single occupancy, and single occupancy 
corresponds in the RFIM to sites with zero spin. Placed 
on a boundary, such sites lower the interaction-energy 
cost of the boundary. We therefore do not expect true 
long-range CDW order in our disordered system. 

Nonetheless, there is a crossover in the vicinity of 
U = AV, where there is a phase transition in the clean 
system. For W < U {V = 0,0.8, and 1.6 in Fig. [T]), 
onsite repulsion remains the dominant interaction, and 
the kinetic-energy driven ZBA found when T^ = per- 



sists. When 4V > U {V = 2.4 and 3.2 in Fig. [T|), the 
nonlocal interaction dominates and the DOS appears to 
be dominated by atomic limit physics. 

The case oi V = 3.2 is particularly simple to under- 
stand from an atomic limit perspective: The DOS has 
a lower band and an upper band separated by a hard 
gap. The location, width and shape of these bands are 
all consistent with the electrons forming a checkerboard 
pattern of alternating empty and doubly occupied sites. 
The lower band corresponds to the removal of particles 
from doubly occupied sites. The nearest neighbors of 
doubly occupied sites are all empty, so the DOS contribu- 
tion of a single site with potential et is Ci + U — fi, with no 
dependence on V . The site potentials e^ are distributed 
between —A/2 and -l-A/2, creating a band in the DOS 
of width A centered on U — ^. The upper band corre- 
sponds to the the addition of particles to empty sites. For 
an empty site there is no onsitc energy cost to adding a 
particle, but the nearest neighbors of empty sites are all 
doubly occupied. In this case, the DOS contribution of a 
single site with potential ci is e^ -|- 8V — fj,. Again, there 
is a distribution of site potentials, resulting in a band of 
width A centered on 8V — jjl. In both cases the DOS 
slants downward toward the Fermi level. This is because 
the ground state for a specific disorder configuration in 
a finite-size system will tend to have sites of especially 
low potential be doubly occupied and sites of especially 
high potential be empty. This means that the lower band 
(coming from doubly occupied sites) will have somewhat 
more contributions from sites with low potentials and 
somewhat fewer contributions from sites with high po- 
tentials. Note that this picture suggests that, although 
true CDW order is not expected, for these parameter val- 
ues the size of the checkerboard domains is larger than 
our system size. 

In this atomic-limit picture, for V — 2.4 the separation 
between the centers of the two bands 8V — U is 11.2, 
slightly less than the disorder strength. The two bands 
have therefore run together, but their slant downward 
toward the Fermi level is still apparent. 



6. Spin and charge correlations 

To further highlight the crossover between y < 1.6 
and V > 2.4, the spin and charge correlations are shown 
in Fig. |4] (a) and (c). For comparison, the spin corre- 
lations for a perfect singly occupied antiferromagnetic 
system are -1-3/4, —1/4 and -1-1/4 for on-site, nearest- 
neighbor, and next-nearest-neighbor respectively. The 
corresponding charge correlations are one for all sepa- 
rations. Meanwhile, the spin correlations for a checker- 
board charge density wave with alternating doubly occu- 
pied and empty sites are all zero, and the charge corre- 
lations are 2, and 2 for on-site, nearest-neighbor and 
next-nearest-neighbor. Essentially these results are con- 
sistent with a crossover from a primarily singly-occupied 
and antiferromagnetically-correlated state at T^ = to 
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FIG. 4: Spin-spin correlation functions as a function of V/t 
(a) at half filling and (b) at quarter filling. Density-density 
correlation functions as a function of V/t (c) at half filling 
and (d) at quarter filling. U = 8t and A = 12i. 



a state with strong charge-density correlations at large 
V. Note that the inflection point in these curves is at 
U ^ W . The simultaneous washing out of the nar- 
row, kinetic-energy-driven anomaly and the suppression 
of nonlocal spin correlations is consistent with the close 
association between these as highlighted in Ref. 1201 . 



7. U dependence at nonzero V 

Fig. [5] shows a series of DOS results at fixed V and 
increasing U . The corresponding points on the phase 
diagram are indicated in the inset. Here again we see a 
crossover from behavior consistent with the atomic-limit 
when U < AV to features uniquely associated with the 
presence of hopping when U > AV. 

For the two lowest U values there is a very broad sup- 
pression of the DOS centered at the Fermi level. This 
is consistent with the picture discussed above of domains 
of checkerboard CDW order. Again, the doubly occupied 
sites in these domains contribute to a plateau in the DOS 
of width A centered at U — fj,, while the singly occupied 
sites generate a corresponding plateau centered at SV—^i. 
Both plateaus slant downwards towards the Fermi level 
due to the tendency for high potential sites to be empty 
and low potential sites to be full in the ground state. In 
the case of f7 = 0, perfect charge ordering would result 
in a very narrow gap, 8V — U — A = 0.8. The presence 
of domain boundaries in a small number of disorder con- 
figurations would fill this in. As U is increased, the gap 
8V — [/ — A closes and in addition the energy cost of 
domain boundaries is lowered. For both these reasons, 
the ZBA is weakened. 

For f/ = 4 and 6, the energy range of the DOS sup- 




FIG. 5: Evolution of the DOS with U/t at half filling for 
A/i = 12 and V/t = 1.6. Inset shows the location of each 
parameter set in the clean phase diagram. 
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FIG. 6: Evolution of the DOS with filling for A/t = 12 
U/t = 8 and V/t = 1.6. 



is related to an important change in the atomic limit 
DOS: For V = Q and in the atomic limit, above quarter 
filling there are doubly occupied sites and below there 
arc none. The Fermi level moves from inside the central 
plateau in Fig. [H^a) to its left edge. As discussed above, 
similar energies of doubly occupied and singly occupied 
states on neighboring sites are central to the emergence 
of the kinetic-energy-driven ZBA. While the presence of 
hopping smooths this transition, by n = 1/4 doping the 
necessary conditions for a linear dependence on t are ab- 
sent. 



pression is sharply reduced. This reduction in the width 
of the ZBA is not associated with hopping t, as it occurs 
in the atomic limit.— An example is seen in Fig. [2][d): 
The curve in Fig. [Ud) shows a ZBA which has roughly 
the same width as those in the U = A and [7 = 6 panels 
of Fig. El 

At U = 8, the onsite interaction is greater than AV 
and the kinetic-energy-driven ZBA unique to strongly 
correlated systems emerges. This ZBA persists through 
U = 12 with a consistent energy scale t, as discussed 
above. 

Finally, a.t U = 16, the Mott gap opens. When V ~ 0, 
the Mott gap opens at C/ ~ A. The addition of non- 
local interactions suppresses single occupancy, delaying 
the formation of the Mott gap and extending the range 
of the strong coupling ZBA. Unlike the Mott gap, the 
strong coupling ZBA is not limited to half filling, and we 
now turn our attention to other dopings. 



B. Away from half filling 

1. Doping dependence 

Figure [H] shows the dependence of the DOS on doping 
for A = 12, C/ = 8 and y = 1.6. The strong similar- 
ities between the curves is consistent with the picture 
that for U > AV there is a kinetic-energy-driven ZBA 
proportional to i. With t and V both held constant, this 
strong coupling ZBA remains unchanged. The deviation, 
especially below the Fermi level, of the n — l/A curve 



2. V dependence at 1/4 filling 
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FIG. 7: Density of states versus frequency at quarter filling 
(n = 1/2) for A = 12i, U = 8t and V values as indicated. 

Figure [7] shows the evolution of the DOS with V at 
quarter filling. Panels (c) and (d) of Fig. |4] show the 
same mean-field rescaling of the DOS discussed for the 
half-filling case in Section IIII A 41 As at half filling, the 
rescaling works well for the V^ = 0.8 case and also at very 
low energies for V = 1.6. 

What is most striking about Fig. [71 however is the 
behavior at large V. The spin and charge correlations 
at quarter filling, shown in Fig. SI (b) and (d), sug- 
gest a crossover similar to the one at half filling, al- 
though much weaker: from AF correlated spins (and no 
CDW correlations) at low V to CDW correlations (and 
no AF spin correlations) at large V. However, whereas 



at half filling (Fig. [T]) these changes in correlations are 
matched by a suppression of the kinetic-energy-driven 
ZBA and the emergence of atomic-limit behavior, at 
quarter filling (Fig. [7|) the form of the ZBA does not 
change when AV > U . Lower filling reduces the compe- 
tition between nearest-neighbor interactions and single 
occupancy. While U reduces double occupancy and V 
suppresses occupation of neighboring sites, both can be 
accommodated by a charge ordered state with alternat- 
ing singly occupied and empty states. Because of the im- 
portance of singly occupied sites to the formation of the 
kinetic-energy-driven ZBA, lower doping may allow it to 
persist to higher V values. It is noticeable, however, that 
persistence of behavior distinct from the atomic limit is 
not confined to the region around the Fermi level. A 
number of other features seen at large V are inconsistent 
with atomic-limit behavior: The onsite charge correla- 
tions increase with V , indicating the presence of dou- 
ble occupancy which is not expected in the atomic limit. 
Also, the DOS remains essentially unchanged by V be- 
low the Fermi level, with no sign of the broad DOS sup- 
pression expected in the atomic limit. Although strong 
interactions reduce the importance of kinetic energy in a 
system, shifting it toward the atomic limit, lower electron 
concentration reduces this eff'ect. 



IV. CONCLUSION 

In conclusion, building on earlier work exploring the 
ZBA in the AHM and its unique strong-coupling fea- 



tures, here we ask how nonlocal interactions influence 
this anomaly. We find that at small values of V and at 
low filling, there is no qualitative change in the anomaly, 
only a gradual renormalization of the hopping amplitude 
which sets its width. At larger values of V and close to 
half filling, however, there is a crossover to DOS features 
which arc independent of hopping. As charge correla- 
tions grow, the energy scale of the kinetic-energy driven 
ZBA gives way to the higher energy scale of the nonlocal 
interaction. The suggestion is that different strongly cor- 
related materials with disorder may display very different 
behaviors depending on where their particular parame- 
ters place them in this phase space. Moreover, it might be 
possible to observe such a crossover in a single material 
by, for example, applying pressure. Increased pressure 
could be expected both to increase the hopping ampli- 
tude and also to reduce the nonlocal interaction through 
increased screening, driving the system from the atomic 
limit toward a regime in which the kinetic-energy driven 
zero-bias anomaly would appear. 
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